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Abstract 

Fuel  control  is  one  of  the  most  pressing  topics  to  achieve  a  self- sustainable  direct  liquid-feed  fuel  cell  system,  such  as  a  direct  methanol  fuel 
cell  (DMFC),  and  enhance  its  overall  efficiency.  In  a  DMFC  system,  sensing  the  methanol  concentration  generally  serves  as  the  basis  of  the  fuel 
control  strategies.  This  paper  proposes  a  three-dimensional  measurement  space  and  constant  concentration  surfaces  (CCS)  to  develop  an  algorithm 
of  estimating  fuel  concentration  in  a  liquid-feed  system,  which  embraces  the  following  merits:  (1)  it  measures  only  three  quantities  or  indices 
that  correlate  with  the  fuel  concentration.  The  indices  can  be  chosen  as  current,  voltage,  temperature,  or  other  quantities  that  are  easily  acquired 
in  an  operating  fuel  cell  system.  The  estimation  can  be  accomplished  without  interrupting  the  operation  of  the  system,  (2)  it  estimates  the  fuel 
concentration  in  a  three-dimensional  measurement  space;  hence  it  is  suitable  for  situations  when  one  or  more  operating  conditions  are  varying, 
(3)  it  can  be  performed  as  a  sensor-less  approach  that  requires  no  additional  methanol  sensors,  thus  consuming  the  less  system  power,  and  (4)  it  is 
particularly  suitable  for  small  and  hand-held  applications. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  suitable  for  small 
or  hand-held  application  systems  like  consumer  electronics 
devices.  In  such  applications,  the  mass,  volume,  and  efficiency 
will  be  the  most  significant  issues  to  be  considered.  However,  the 
methanol  crossover  phenomenon  in  a  DMFC  has  been  a  crucial 
limitation  that  may  greatly  reduce  the  efficiency  of  the  system 
when  methanol  fuel  of  higher  concentration  is  fed  into  the  fuel 
cell  units.  On  the  other  hand,  to  achieve  a  specified  total  energy 
delivery,  introducing  methanol  fuel  of  lower  concentration  will 
lead  to  a  lower  energy  density;  and  the  system  would  be  heav¬ 
ier  and  larger.  To  take  the  efficiency  and  energy  density  into 
account  simultaneously,  it  is  desirable  to  replenish  a  fuel  car¬ 
tridge  with  concentrated  methanol  fuel  to  increase  the  energy 
density  of  the  system.  Then  the  fuel  is  pre-diluted  inside  the 
fuel  cell  system  to  an  adequate  concentration  level  before  being 
fed  onto  the  fuel  cell  units  (i.e.  membrane  electrode  assemblies, 
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ME  As),  thus  preventing  excessive  fuel  crossover.  Furthermore, 
fuel  concentration  is  one  of  the  most  significant  factors  that 
dominate  the  behavior  of  a  fuel  cell.  Hence  it  should  be  well 
controlled  to  achieve  optimal  performance  in  a  fuel  cell  system. 
As  a  result,  fuel  concentration  control  is  essential  to  meet  the 
requirements.  Techniques  of  sensing  methanol  concentration  are 
therefore  important  to  serve  as  the  basis  of  fuel  control  strategies. 

There  have  been  various  approaches  to  acquire  the  methanol 
concentration  in  a  DMFC  system.  The  approaches  can  be  divided 
into  two  major  categories.  The  first  one  correlates  the  con¬ 
centration  of  the  methanol  aqueous  solution  with  its  physical 
properties,  such  as  density,  dielectric  constant,  ultrasound  veloc¬ 
ity,  infrared  absorptivity,  optical  behavior,  transmittance,  etc. 
[1-5].  However,  this  kind  of  approach  requires  more  complex 
and  delicate  instruments  that  may  not  be  suitable  for  small  and 
hand-held  applications  under  cost,  mass,  and  volume  constrain. 
The  second  category  of  methanol  concentration  sensing  schemes 
is  based  on  the  electrochemical  characteristics  of  the  fuel  cells 
themselves.  They  can  further  be  divided  into  two  types:  galvanic 
type  and  electrolytic  type.  The  former  correlates  the  polarization 
curve  [6],  short  circuit  current  [7],  open  circuit  voltage  [8],  and 
even  indirect  measurement  data  [9,10]  with  the  methanol  con- 


Y.-J.  Chiu,  H.-C.  Lien  /  Journal  of  Power  Sources  159  (2006)  1162-1168 


1163 


centration  to  be  determined.  These  sensing  schemes  adopt  the 
inherent  behavior  of  a  fuel  cell  and  hence  are  simple,  straight¬ 
forward,  and  even  sensor-less.  However,  they  are  very  sensitive 
to  the  operational  conditions  such  as  temperature,  voltage  (or 
current  density),  and  flow  rate  of  the  fuel  (or  air).  Additional 
compensation  is  needed  for  situations  that  are  varying.  The  elec¬ 
trolytic  type  sensing  schemes  measure  the  limiting  current  or  the 
half-cell  potential  to  establish  the  relationship  with  methanol 
concentration  [11-15].  In  these  designs,  a  substantial  sensor  is 
required  while  the  reactions  wherein  are  not  spontaneous.  Addi¬ 
tional  external  power  input  is  required  to  drive  the  oxidation 
and  reduction  reactions  on  the  anode  and  cathode  sides,  respec¬ 
tively.  Temperature  still  plays  a  significant  role  and  influences 
the  measurement  results. 

The  present  paper  proposes  an  estimation  algorithm  to  deter¬ 
mine  methanol  concentration  in  a  small  DMFC  system.  The 
method  is  based  on  a  three-dimensional  measurement  space 
with  the  measurement  dimensions  being  current,  voltage,  and 
temperature.  Whenever  other  measurement  quantities  or  indices 
can  be  correlated  with  the  methanol  concentration,  they  could 
also  serve  as  the  three  measurement  dimensions  while  the  pro¬ 
posed  algorithm  is  still  valid.  Since  the  method  is  performed 
in  a  three-dimensional  measurement  space,  it  is  suitable  for  sit¬ 
uations  where  one  or  more  operating  factors,  like  temperature 
or  current  density,  are  varying.  Moreover,  it  can  be  designed 
as  a  sensor-less  approach  that  consumes  the  less  system  power 
while  minimizing  the  cost,  mass,  and  volume  of  the  system. 
Such  a  design  is  particularly  beneficial  for  small  and  hand-held 
applications. 

2.  The  estimation  algorithm 

The  present  paper  proposes  an  Interpolation  algorithm  based 
on  Constant  Concentration  Surfaces  (ICCS  algorithm)  to  esti¬ 
mate  the  fuel  concentration  in  a  liquid-feed  fuel  cell  system. 
Its  primary  concept  accompanying  the  solution  schemes  will  be 
illustrated  in  this  section. 


2.7.  The  notion  of  ICCS  algorithm 


The  polarization  curve  (I-V  curve)  reveals  the  most  impor¬ 
tant  characteristic  of  a  fuel  cell,  which  depicts  the  relation  of 
voltage  to  current  density.  The  relation  can  be  briefly  expressed 
as  follows  [16]: 
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where  E  denotes  the  reversible  open  circuit  voltage,  while  V  and 
i  are  voltage  and  current  density,  respectively.  For  a  given  direct 
methanol  fuel  cell,  the  values  of  the  parameters  A,  B ,  io,  in,  // , 
and  r  in  Eq.  (1)  closely  depend  on  the  operational  at  conditions 
including  temperature,  methanol  concentration  supplied  to  the 
anode,  flow  rates  of  the  delivered  fuel  and  air,  and  pressure  built 
in  the  anode  and  cathode  chambers.  The  delivery  of  fuel  and 
air  not  only  provides  the  reactants  for  the  electrochemical  reac¬ 
tions,  but  also  removes  the  products  and  heat  generated  by  the 
reactions.  Accordingly,  the  effect  of  flow  rates  of  fuel  and  air 


on  Eq.  (1)  can  be  indirectly  translated  into  effects  on  the  tem¬ 
perature  and  fuel  concentration.  Moreover,  the  pressure  effect 
can  influence  the  mass  transport  or  concentration.  Generally, 
the  pressure  issue  is  omitted  in  a  smaller- scaled  DMFC  system 
since  only  a  limited  energy  budget  is  left  for  the  balance  of  plant 
(BOP),  which  cannot  offer  enough  power  to  build  a  much  higher 
pressure  for  fuel  and  air  delivery.  As  a  result,  whenever  the  flow 
rates  of  fuel  and  air  are  specified,  the  temperature  and  methanol 
concentration  will  be  the  most  significant  condition  variables 
to  be  controlled.  Then  for  a  given  DMFC,  the  behavior  of  the 
polarization  curve  in  Eq.  (1)  can  be  written  as 

/(/,  F,  C,  T)  =  0  (2) 

where  /,  V,  C  and  T  denote  current,  voltage,  fuel  concentration, 
and  operating  temperature  in  order.  One  can  observe  that,  in 
Eq.  (2),  three  of  the  four  variables  are  independent  while  the 
rest  should  be  compliant  with  the  equation.  It  represents  the 
most  important  relationship  between  the  material  properties  and 
system  operating  variables. 

Generally  speaking,  the  aforementioned  operational  vari¬ 
ables  should  be  confined  within  certain  regions,  e.g.  20-70  °C 
and  0.5-2.0M  for  temperature  and  fuel  concentration,  respec¬ 
tively.  Under  such  situations,  the  I-V  curve  of  a  DMFC  will 
positively  shift  toward  the  higher  performance  region  (higher 
voltage  with  fixed  current  or  higher  current  with  fixed  volt¬ 
age)  whenever  the  operation  temperature  increases  (adopting 
constant  fuel  concentration)  or  the  fuel  concentration  increases 
(adopting  constant  operation  temperature).  Although  the  afore¬ 
mentioned  qualitative  tendency  of  the  DMFC  properties  has 
been  well  known,  it  is  still  a  complicated  task  to  quantitatively 
determine  the  relationship  described  in  Eq.  (2). 

Since  three  of  the  four  variables  in  Eq.  (2)  are  independent, 
this  paper  introduces  a  three-dimensional  measurement  space 
m3  spanned  by  three  measurement  axes  as  its  basis  to  estimate 
fuel  concentration.  In  the  case  of  Eq.  (2),  if  the  fuel  concentra¬ 
tion  level  is  specified,  say  Co,  then  a  series  of  I-V  curves  can  be 
acquired  corresponding  to  various  temperature  levels.  The  I-V 
curves  are  depicted  in  a  three-dimensional  measurement  space 
91 3  with  its  axes  being  current,  voltage,  and  temperature.  More¬ 
over,  a  surface  S  spanned  by  these  I-V  curves  can  be  obtained, 
as  is  shown  in  Fig.  1.  Since  all  these  I-V  curves  correspond  to 
the  same  fuel  concentration  level  Co  and  they  will  shift  posi¬ 
tively  toward  the  higher  performance  region  accompanying  an 
increase  of  temperature,  any  other  point  on  this  surface  S  will 
also  be  correlated  to  Co.  Accordingly,  such  a  surface  is  defined 
as  a  Constant  Concentration  Surface  (CCS). 

For  an  unknown  fuel  concentration  level  Cu,  one  can  take 
an  in  situ  measurement  of  current,  voltage,  and  temperature  of 
a  fuel  cell  system  then  have  a  point  P+  (or  P“)  located  in  the 
measurement  space  9t3,  as  was  shown  in  Fig.  2.  By  judging  how 
far  the  point  P+  (or  P_)  is  away  from  the  surface  S,  it  is  inferred 
how  large  the  concentration  difference  between  Cu  and  Co  is. 
Whenever  two  or  more  CCSs  are  pre-established  with  given 
fuel  concentration  levels,  an  Interpolation  algorithm  based  on 
Constant  Concentration  Surfaces  (ICCS)  can  then  be  adopted 
for  estimating  the  unknown  concentration  level  Cu,  as  will  be 
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2.2.1.  Establishment  of  C CSs 

First,  a  number  of  CCSs  are  pre-established  corresponding  to 
specified  fuel  concentration  levels.  Each  CCS  is  approximated 
by  a  cubic  polynomial  equation.  Therefore,  only  10  coeffi¬ 
cients  should  be  stored  in  the  energy  management  system  (EMS) 
instead  of  the  whole  experimental  database,  thus  greatly  enhanc¬ 
ing  its  efficiency. 

By  confining  a  specified  fuel  concentration  level  Q,  a  set  of 
I-V  curves  can  be  obtained  experimentally  under  various  tem¬ 
perature  levels  and  a  CCS  S^  is  generated  in  9f3.  The  surface  S^ 
is  defined  as  follows 


Fitted  CCS  of  5v%  MeOH 


voltage  (V)  0  6  30  temperature  (°C) 


Fig.  1.  A  constant  concentration  surface  in  a  three-dimensional  measurement 
space  . 

illustrated  later.  Consequently,  the  present  paper  depicts  Eq.  (2) 
by  a  graphical  manner,  meanwhile,  gives  a  well  estimation  of 
fuel  concentration  in  a  liquid-feed  fuel  cell  system. 

2.2.  The  solution  scheme  of  ICCS  algorithm 

The  solution  scheme  of  the  ICCS  algorithm  can  be  divided 
into  two  stages,  that  is,  establishment  of  CCSs  and  the  interpo¬ 
lation  procedure. 


fk(I,  V. ;  C,  T)\c=ck  =  gM  v,T)  =  0  (3) 

For  any  combination  of  (/,  V,  T)  that  satisfies  Eq.  (3),  it  repre¬ 
sents  a  point  on  the  surface  Sk  corresponding  to  the  same  fuel 
concentration  level  Q.  In  practice,  let  function  gk  in  Eq.  (3)  be 
approximated  by  a  cubic  polynomial  as 

gk(f  V,  T)  =  I  —  (< a\k  +  a2kT  +  a^k V  +  a^T2  +  a^kTV 

+  a6k  V2  ~\~alk  T2  -\~a%k  T2  V  +  agkTV2  +  <3io  &V3) 

=  0  (4) 

where  a\k,  a2k ,  . . <3io£  are  constant  coefficients  to  be  deter¬ 
mined.  One  can  take  measurement  (//,  Vi,  Tf)  for  i=  1,  2,  . . .,  m 
under  various  temperature  and  voltage  (or  current)  levels  and 
substitute  ah  the  values  into  Eq.  (4)  to  have  a  set  of  polynomial 
equations  consisting  of  unknown  coefficients  a\k,  a2k ,  . . 
<310  k' 
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Sampled  Points  and  CCS 


Fig.  2.  A  pre-established  CCS  and  measurement  points  in  9f3. 


Then  the  well  known  Least-Squares  Method  can  be  adopted  for 
solving  ak  =  [ak,  a\k,  a2k, .  ■ awk]T  in  Eq.  (5): 

ak  =  (MTM)_1MTb  (6) 

A11  the  data  acquisition  and  relevant  computations  can  be  per¬ 
formed  beforehand,  while  only  the  values  of  ak  are  necessary  in 
the  energy  management  system  of  a  fuel  cell  system  to  represent 
a  pre-established  CCS  of  concentration  Q. 

2.2.2.  Interpolation  procedure  based  on  CCS 

Let  an  arbitrarily  point  Pu(/U,  Vu,  Tu)  in  9t3  correspond  to  an 
unknown  fuel  concentration  level  Cu,  where  /u,  Vu,  Tu  are  the 
measured  current,  voltage  and  temperature,  in  that  order.  More¬ 
over,  a  number  of  CCSs  S k  corresponding  to  the  given  methanol 
concentration  levels  Ck,  for  k=  1,  2,  . . .,  n,  are  pre-established. 
The  following  illustration  will  transfer  the  relationship  between 
Pu(Ai>  Vu,  Tu)  and  S^  in  9f3  into  a  current-concentration  relation 
to  have  an  interpolation  formula. 

Whenever  a  number  of  CCSs  S k  and  a  measured  point  Pu(Ai, 
Vu,  Tu)  are  depicted  in  91 3 ,  Pu  will  have  n  projection  points  Q k 
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on  S&  for  k=  1,  2,  . . n  along  the  current  axis.  The  projection 
point  Qk  has  a  value  of  current  4  that  can  be  interpreted  as 

4  =  <2\ k  +  CllkJw  +  a?>k  Vu  +  (24  kT^  +  (25kTuVu  +  CLfrkVu 

+  (27  k  4i  +  a&k  Vu  +  agk  TUV 2  +  (2 10  k  V y  (7) 


where  Vu  and  Tu  are  the  voltage  and  temperature  of  the  measured 
point  Pu,  while  a\k ,  a^k ,  flioit  are  the  coefficients  of  the  pre- 
established  CCS  of  concentration  Q  for  k-  1,  2,  AT, 

Then  the  unknown  fuel  concentration  Cu  can  be  determined 
by  the  Lagrange  interpolation  formula  as  follows  [17]: 


n 


cu  =  E 

k=  1 


where  7U  is  the  current  value  of  the  measured  point  Pu,  Ck  denotes 
the  fuel  concentration  of  the  pre-established  SSCs,  while  Ij  and 
4  is  determined  by  Eq.  (7). 

Accordingly,  the  ICCS  algorithm  has  an  in  situ  and  unique 
estimation  of  the  fuel  concentration  of  a  liquid-feed  fuel  cell  sys¬ 
tem,  merely  by  measuring  the  current,  voltage  and  temperature 
that  are  all  easily  acquired  throughout  the  operating  region.  One 
remark  should  be  made  that  whenever  other  measurement  quan¬ 
tities  or  indices  can  be  correlated  with  the  methanol  concentra¬ 
tion,  they  could  also  serve  as  the  three  measurement  dimensions 
instead  of  current,  voltage  and  temperature,  and  the  proposed 
algorithm  is  still  valid. 


3.  Experiments 

The  aforementioned  ICCS  algorithm  is  verified  by  experi¬ 
mental  results  based  on  standard  testing  procedures  of  a  unit 
fuel  cell.  A  standard  testing  module  of  a  unit  DMFC  (with 
50  mm  x  50  mm  MEA,  DuPont  Nation®  117)  is  adopted  on 
a  testing  system  (890B  model  of  Scribner  Associates  Inc.)  in 
the  following  experiments.  The  channel  type  of  the  graphite 
endplates  is  serpentine  and  2  mm  in  width  and  depth.  Four  tem¬ 
perature  levels  and  six  concentration  levels  are  considered  in 
establishing  CCSs,  as  are  shown  in  Table  1.  The  methanol  fuel 


Table  1 

A  brief  list  of  the  experimental  conditions 


Items 

Descriptions 

1.  MEA 

50  mm  x  50  mm  (DuPont  Nation®  117) 

2.  Graphite  channel  type 

Serpentine,  2  mm  in  width  and  depth 

3.  Temperature  levels  (°C) 

30,  40,  50,  60 

4.  Fuel  concentration  levels  for 

3,  4,  5,  6,  8,  10 

pre-established  CCSs  (vol.%) 

5.  Fuel  concentration  levels  for 

3.5,  4.5,  5.5 

verification  points  (vol.%) 

6.  Fixed  air  flow  rate  (ml  min  1 ) 

500 

7.  Fixed  fuel  flow  rate 

1 1 .2  ml  min~ 1  (2 A  at  4  x  stoichiometric 
rate  in  3  vol.%  fuel) 

8.  Electrical  load 

Scanned  under  constant  voltage  mode 
from  0.5  to  0. 1  V  with  an  increment  of 

-0.05  V 

is  pre-heated  to  the  desired  temperature  then  delivered  by  a  liquid 
pump  and  accompanying  flow  rate  control.  Moreover,  a  thermo¬ 
stat  is  attached  to  the  testing  module  to  control  the  operating 
temperature  during  the  experiments.  In  order  to  keep  the  fuel 
concentration  constant  in  each  level,  the  fuel  exhausted  from 
the  outlet  of  the  testing  module  will  not  circulate  back  to  the 
fuel  reservoir.  Furthermore,  the  flow  rates  of  the  methanol  and 
air  are  kept  constant  to  exclude  the  influence  of  heat  exchange 
deviation.  That  would  be  500  ml  min  - 1  for  air  and  1 1 .2  ml  min  - 1 
for  the  methanol  solution  (2 A  at  4x  stoichiometric  rate  in  3% 
methanol  fuel).  For  a  given  temperature  and  concentration,  the 
electrical  load  is  controlled  under  a  constant  voltage  mode  from 
0.5  to  0. 1  V  with  increments  of  —0.05  V.  Meanwhile  an  I-t  chart 
(current  versus  time)  is  recorded.  There  will  be  corresponding 
rising  steps  in  the  I-t  charts  whenever  the  voltage  is  downward 
regulated.  The  stabilized  current  values  in  each  step  are  acquired 
for  each  voltage  scan  to  have  an  I-V  curve.  For  each  concentra¬ 
tion  level,  there  will  be  four  I-V  curves  of  different  temperature 
levels  and  36  measured  points  of  (4,  Vi,  Tf)  in  total,  thus  estab¬ 
lishing  the  CCS  by  using  Eqs.  (5)  and  (6). 

The  following  experimental  evaluation  consists  of  three  parts: 
(1)  concept  validation  of  CCS,  (2)  establishment  of  CCSs  and 
feasible  measurement  regions,  and  (3)  evaluation  of  ICCS  algo¬ 
rithm. 

4.  Results  and  discussions 

4. 1.  Concept  validation  of  CCS 

Before  proceeding  with  further  evaluation,  the  concept  of 
CCS  is  first  validated.  Four  I-V  curves  of  different  temperature 
levels  listed  in  Table  1  are  acquired  based  on  5  vol.%  concen¬ 
tration.  Then  the  constant  coefficients,  denoted  by  as,  can  be 
solved  by  Eqs.  (5)  and  (6),  thus  establishing  a  CCS  named  S5. 
Afterwards  16  testing  points  P )(//,  V/,  7)),  for  j=  1,  2,  . . .,  16, 
randomly  spread  over  the  temperature  interval  of  30-55  °C  and 
the  voltage  interval  of  0. 1-0.4  V  are  measured,  again  with  the 
methanol  fuel  at  5  vol.%  concentration.  By  substituting  Vu  =  Vj, 
Tu  -Tj  and  as  into  Eq.  (7),  one  can  obtain  current  value  Ij  of 
the  projection  point  Qj  on  S 5.  If  Ij  =  Ij,  it  indicates  that  P )  is 
on  the  CCS  S5  that  satisfies  the  definition  of  CCS.  Fig.  3  shows 
the  CCS  S5  and  the  16  testing  points,  where  symbols  ‘O’  and 
‘+’  represent  P y  and  Qj,  respectively.  The  deviations  between  Py 
and  Qj,  defined  as  (Ij  ~  7j)/7j  x  100%,  are  depicted  in  Fig.  4. 
It  turns  out  that  all  the  deviations  of  the  testing  points  are  less 
than  5%,  thus  sufficiently  validating  the  notion  of  CCS. 

4.2.  Establishment  of  CCSs  and  feasible  measurement 
regions 

Following  a  similar  procedure  that  creates  S5,  the  other  five 
CCSs  corresponding  to  the  concentration  levels  listed  in  Table  1 
are  established.  The  experimental  results  show  that  CCSs  may 
intersect  with  each  other  due  to  the  inherent  characteristics  of 
a  DMFC.  Such  a  situation  will  lead  to  ambiguities  in  the  ICCS 
algorithm,  hence  the  feasible  regions  of  concentration  and  volt¬ 
age  should  be  well  defined. 
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Sampled  Points  and  CCS  with  5v%  MeOH 


voltage  (V)  0.6  30  temperature  (°C) 


Fig.  3.  Pre-established  CCS  S5,  individually  sampled  points  P;,  and  the  projec¬ 
tion  points  Qj  with  the  same  concentration  level  (5  vol.%). 


I-V  Curves  of  Different  Concentration  Levels  (at  30°C) 
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Fig.  5.  The  I-V  curves  of  different  concentration  levels  (at  30  °C). 


First,  the  I-V  curves  with  different  concentration  levels  are 
depicted  in  Fig.  5  (at  30  °C).  It  shows  that  all  the  curves  inter¬ 
sect  with  each  other  near  the  region  of  higher  voltage  (or  lower 
current  density).  Meanwhile  in  the  so-called  mass  transport  or 
concentration  losses  region  (generally  at  lower  voltages),  the  I-V 
curves  spread  out  by  the  increase  of  concentration  in  an  ordered 
manner.  As  a  result,  it  is  less  to  employ  the  estimation  algorithm 
in  the  lower  voltage  region  less  than  0.35  V  to  avoid  ambiguities 
and  achieve  increased  accuracy  in  the  results.  However,  owing  to 
efficiency  consideration,  it  is  not  advisable  to  operate  a  DMFC 
below  the  lower  voltage  region.  Consequently,  for  a  larger-scale 
DMFC  system,  it  is  suggested  that  a  sensor  unit  is  included 
separate  from  the  fuel  cell  stack  that  generates  electricity.  The 
two  parts  can  operate  independent  at  different  voltage  intervals. 
As  for  a  smaller-scale  DMFC  system,  to  simplify  the  system 
architecture,  a  sensor-less  approach  can  be  used.  The  fuel  cell 
operates  at  a  higher  voltage  in  the  normal  state.  Whenever  the 


Deviation  of  the  Sampled  Points  from  the  Fitted  CCS  (5v%  MeOH) 
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Fig.  4.  Relative  error  (deviation  between  P )  and  Q;  )  in  current  direction  regarding 
the  CCS  S5. 


concentration  information  is  needed,  the  energy  management 
system  lowers  the  voltage  of  the  fuel  cell  to  acquire  the  neces¬ 
sary  measurement  data,  and  then  returns  to  the  normal  state.  In 
the  present  example,  the  time  required  for  taking  a  stabilized 
measurement  point  Pu(/u,  Vu,  Tu)  is  about  2-6  s. 

Fig.  5  also  reveals  that  when  the  concentration  is  within  the 
region  of  3-8  vol.%,  the  I-V  curves  spread  out  with  increase  of 
concentration  in  an  ordered  manner.  However,  when  the  concen¬ 
tration  is  further  increased  to  10  vol.%,  the  curve  intersects  the 
others  and  shows  a  tendency  to  move  ‘backward’ .  This  is  mainly 
due  to  the  methanol  crossover  phenomenon  and  can  become 
more  significant  at  higher  temperatures.  Accordingly,  the  fea¬ 
sible  region  of  the  measured  concentration  should  be  confined 
within  a  certain  range.  In  the  present  example,  it  should  not 
exceed  8  vol.%  (2.0  M).  Fortunately,  in  a  practical  DMFC  sys¬ 
tem,  methanol  crossover  should  be  properly  suppressed  to  ensure 
adequate  fuel  utilization  efficiency  and  Faradaic  efficiency.  The 
literature  has  indicated  that  the  suggested  fuel  concentration 
would  be  less  than  1.5-2.0M  (6-8  vol.%)  [18,19].  In  the  fol¬ 
lowing  illustration,  the  methanol  concentration  will  be  confined 
to  3-6  vol.%  while  the  corresponding  CCSs  are  shown  in  Fig.  6. 

The  above  limitations  are  derived  by  the  inherent  characteris¬ 
tics  of  a  DMFC.  All  the  electrochemical  type  methanol  sensors 
will  more  or  less  suffer  from  similar  situations.  The  feasible 
regions  depend  on  the  MEA  design  of  a  fuel  cell  system.  For 
example,  a  membrane  coated  with  a  thin  palladium  (Pd)  film 
will  dramatically  reduce  the  methanol  crossover  phenomenon 
[5],  thus  extending  the  range  of  the  measured  concentration.  In 
such  a  situation,  the  present  ICCS  algorithm  will  also  possess 
a  larger  measurement  domain.  Furthermore,  if  other  measure¬ 
ment  indices  are  chosen  as  the  basis  of  9i3,  the  ranges  of  the 
measurement  practice  may  then  have  corresponding  changes. 

4.3.  Evaluation  of  ICCS  algorithm 

To  evaluate  the  proposed  ICCS  algorithm,  three  additional 
concentration  levels  (shown  in  Table  1)  are  introduced  to  serve 
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Fitted  CCSs  in  the  Feasible  Region 


temperature  (°C) 


0  15  0.2  0.25  0.3 

voltage  (V) 


0.35 


Estimation  Error  of  the  ICCS  Method 


Fig.  6.  The  CCSs  in  the  feasible  region  that  guarantees  the  ICCS  algorithm  to 
have  a  unique  and  robust  estimation. 


Fig.  8.  The  evaluation  of  the  proposed  ICCS  algorithm  9i3. 


as  the  unknown  concentration  to  be  estimated.  Each  of  the  three 
levels  is  assigned  12  measured  points  randomly  in  the  temper¬ 
ature  range  from  30  to  55  °C  with  a  5  °C  increment,  and  the 
voltage  range  from  0.1  to  0.3  V  with  a  0.05  V  increment.  The 
distribution  of  the  36  testing  points  in  the  three-dimensional 
measurement  space  91 3  is  shown  in  Fig.  7.  Since  there  are  four 
pre-established  CCSs  (in  Fig.  6),  the  ICCS  algorithm  will  have 
a  3rd  ordered  estimation  when  adopting  Eqs.  (7)  and  (8).  The 
results  are  shown  in  Fig.  8.  It  indicates  that  the  estimation  is 
quite  robust  since  all  the  estimated  errors  have  distributed  over  a 
uniform  range  of  ±0.3  vol.%  (about  ±0.075  M).  This  is  believed 
to  have  adequate  accuracy  for  a  real  application. 

In  the  aforementioned  illustration,  the  proposed  approach 
provides  good  results  under  the  situations  that  both  the  tem¬ 
perature  and  voltage  are  varying.  Hence  it  is  more  suitable  for 
in  situ  measurement  where  the  impact  of  the  surroundings  may 


Sampled  Points  in  the  3D-Measurement  Space 


Fig.  7.  The  distribution  of  the  testing  points  in  the  three-dimensional  measure¬ 
ment  space  9± . 


influence  the  operation  conditions  of  the  system.  Accordingly, 
the  feasibility  and  accuracy  of  the  proposed  strategy  are  verified. 

5.  Conclusions 

A  strategy  for  estimating  the  methanol  concentration  in  a 
DMFC  is  developed.  The  approach  gives  a  clear  measurement 
of  electrochemical  behavior  of  a  DMFC  in  Eq.  (2)  by  intro¬ 
ducing  the  notion  of  CCS  in  a  graphical  manner.  It  is  also  an 
efficient  method  to  record  a  set  of  coefficients  depict  a  CCS  in 
the  energy  management  system  of  a  fuel  cell  system,  instead  of 
the  whole  experimental  database.  In  practice,  the  in  situ  esti¬ 
mation  performs  only  Eqs.  (7)  and  (8)  instead  of  looking  table 
schemes  in  a  multi-dimensional  space.  Therefore,  it  is  suitable 
for  dynamic  operating  conditions,  taking  surrounding  influences 
such  as  temperature  into  account. 

The  primary  concept  of  three-dimensional  measurement 
space  and  a  constant  concentration  surface  have  been  experi¬ 
mentally  verified.  The  proposed  algorithm  is  shown  to  possess 
good  results  over  the  range  of  practical  operating  conditions.  It 
is  also  suitable  for  other  kinds  of  liquid-feed  fuel  cells  where 
the  fuel  concentration  should  be  measured.  Since  the  method 
can  be  performed  in  a  sensor-less  approach  that  requires  none 
of  additional  methanol  sensors,  it  is  particularly  useful  for  small 
hand-held  applications. 
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